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Abstract 

Numerical simulations of the scattering of a linear plane wave incoming onto a nonlin- 
ear medium (sine-Gordon) reveals that: i) nonlinearity allows energy transmission in the 
forbidden band, ii) this nonlinear transmission occurs beyond an energy threshold of the 
incoming wave, iii) the process begins (at the threshold) with large amplitude breathers, 
and then energy is generically transmitted both by kink-antikink pairs and breathers. 

Introduction From the Fermi- Pasta- Ulam recursion phenomenon 0, nonlinear lattice and 
their numerical simulations have played a key role in the development of nonlinear studies. 
Since then a whole branch of theoretical physics has evolved with the birth of the concepts 
of soliton H and integrability || , for continuous, discrete and higher dimensional systems M . 
These discoveries have had a enormous impact on the understanding of nonlinear physics ||. 



Q\ ', One of the fundamental studied problem is the energy localization (and transport) by gener- 
O ation (and propagation) of nonlinear coherent structures (solitons, breathers, nonlinear modes) 
^ 'by many different processes such as external forcing, noise, local potential, inhomogeneous forc- 
O ' 



O 



ing, anharmonic perturbation, modulated external forcing etc... [H]. We note that nowadays 
this problem is more largely studied in the context of discrete systems where intrinsic nonlinear 
modes (breathers for short) have been shown to play the main part in storage and transport of 
energy [0. 

Another basic question is the scattering of waves incoming onto a nonlinear medium. Two 
'classes of problems related to wave scattering have been considered in detail: the resonant 
^ .wave coupling which includes self-induced transparency (SIT) [|J, stimulated Raman scattering 



two-photon propagation [11], etc., and the self modulation resulting from Kerr effect in 
optical media [12|] ||13|| . Within these two classes of problems, solitons have been shown to play 
the central role in propagating energy in regions of parameter (frequency, energy,...) usually 
forbidden by the linear theory. 

The resonant wave coupling describes the process of interaction of an electromagnetic radi- 
ation with a two- level medium through Brillouin selection rules (the nonlinearity is extrinsic). 
For SIT the linear theory would conclude with resonant absorption of radiation while nonlinear- 
ity allows for perfect transparency. The model equations read as coupled Maxwell-Bloch system 
and the scattering of waves results in a boundary-value problem for the reduced system. The 



point there is that this boundary-value problem is integrable - see also [14] - and it has been 
shown that the soliton (actually a kink for the phase) is the means which propagates energy 
through the medium without loss or distortion, provided the input pulse obeys the area theorem 
||. Hence in SIT, the soliton is created by a sufficiently energetic boundary condition. In the 
inverse spectral transform scheme, this boundary value problem actually maps to a Cauchy 
problem and the area theorem simply states that the initial datum must produces at least one 



discrete eigenvalue || |T~3| . 
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In the case of self modulation, the nonlinearity finds its origin in the medium itself (intrinsic 
Kerr effect). A process similar to SIT (corresponding to a Cauchy problem in the 1ST scheme), 
but with different origin, occurs in fiber guides in nonlinear regime (and with anomalous dis- 
persion) where a sufficiently energetic input light pulse evolves as a soliton of the nonlinear 
Schrodinger equation (NLS) along the fiber [|TJ] |I2| . Here and in SIT, the energy flows with 
solitons created by a localized high energy pulse. 

A different soliton effect arises in the study of the propagation of light in a medium with 
forbidden bands (photonic band gaps) for which the concept of gap soliton have been introduced 
||16|| . There it has been shown that the nonlinearity in a dielectric periodic thin-film stack 
(medium with spatial periodic variation of the refraction index) allows for transmission of 
radiation inside the gap by means of an envelope soliton. In short, the presence of an envelope 
soliton in the medium shifts locally the gap and allows for transparency [|17]]. The model 
for the slowly varying envelope is the NLS equation fl8f . For that reason (transmission of 
radiation in the gap), and the related switching properties (tranmission with soliton - reflection 



without), photonic band gap media are the subject of active research, see e.g. Notice 
that the creation of the envelope soliton in those media is an open question. Indeed, while 
the nonlinearity, via NLS, offers possible existence of a soliton, the mechanism of its generation 
itself is still to be understood. Such a mechanism has been for instance analytically described in 
the case of stimulated Raman scattering in ]10[ with the tool of the inverse spectral transform 
applied to boundary value problems (another example can be found in [pCfl ). 



Then the mechanism of soliton creation by radiation is of fundamental interest and we 
consider here the (nonlinear) scattering of waves from a quite elementary point of view. Namely 
we explore the possibility for a nonlinear medium to build up nonlinear coherent structures 
under "exposition" to input radiation. The problem we are interested in, very simple in linear 
systems, concerns the reflection and transmission of a plane wave in an inhomogeneous medium 
with band pass type dispersion relation: a plane wave created in a medium 1 is sent onto a 
medium 2 with a different dispersion relation. In particular total reflection occurs when the 
incoming wave of medium 1 falls onto medium 2 with a frequency inside a frobidden band. We 
consider then this generic question in the case when the medium 2 is intrinsically nonlinear and 
modelized by an integrable equation (sine- Gordon). 



Statement of the problem. We consider then two different media defined by the wave 
equations 

x < : u tt - c\u xx + u]\u = , (1) 

x > : u t t — c\u xx + uj\ sin u = , (2) 

with continuity conditions in x = (inferred here from the natural discrete version of the above 
system, see later). The scattering of a wave of given frequency Q and normalized amplitude 
incoming from medium 1 to medium 2, when Q lies in the phonon gap of the nonlinear medium 
(fl < U2), is described in the linear limit by the general solution (the normalization is chosen 
such as to have an incident wave Asin(Kx — Qt)) 

x < : u(x, t) = ^-.e- int [Ae iKx + Be~ iKx ] + c.c. , (3) 

x > : u(x, t) = ^e~ iQt [Ce~ KX + De RX ] + c.c. , (4) 

with K and k given from Q by 

n 2 = uj\ + K 2 c\ = w| - k 2 c 2 2 . (5) 



While the linear theory produces an evanescent wave (D = 0), the nonlinearity does not 
prevent in principle nonlinear superpositionf] of exponentials in (Q), possibly allowing energy 
transmission to the nonlinear medium by means of a soliton. We address then the question 
wether or not solitons build up under scattering and allow nonlinear transmission of radiation 
within the gap. 

The present study is fully based on numerical simulations and, before going further, it is 
worth stressing that we will discover that indeed solitons (actually kink-antikink pairs and 
breathers) are created to allow for energy transmission. Moreover, this process is in essence 
different either from self-induced transparency in two-level media, or else from gap solitons in 
photonic band gap media, because reductive perturbation method (or multi-scale expansions) 
do not apply: the incoming radiation creates large amplitude nonlinear objects. 



Method. The system (fj])(0) is first maped to the following discrete analog where the first 
nonlinear pendulum U\ is coupled to the last linear pendulum uq by the spring c\ of the linear 
chain: 
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u n - c\(u n+ i - 2u n + u n ^i)/h 2 + io\u n = 



cl(u 2 — Mi) /h 2 + c 2 (ui — Mo) /h 2 + u)\ sin u\ = , 

cl(u n+ i — 2u n + M n _i)//i 2 + Uj\ sinM n = . (6) 

This maps to (P(||) in the continuous limit h — > and has natural unconstrained continuity 
conditions in n = 1. 

In order to generate in the linear medium 1 an incoming plane wave we follow the work |2T 
by making the change of variable from u n (t) to f n (t) : 

Un (t) = f n (t) + AP{n) sm{Knh - ttt) , (7) 

where Q is given from K by (^), where A is the amplitude of the incoming wave and where the 
fuction P(n) (a polynomial) is chosen such as to vary from the value 1 for large negative n to 
the value far before the nonlinear medium 2. To set things clear, we have used a medium 1 
made of a (linear) chain of 10500 oscillators and the plynomial P(n) varies from the value 1 at 
n = —10500 to zero around n = —3000. The medium 2 (nonlinear chain) is constituted of 4500 
pendula. Finally, to simulate the infinite line, strong damping is included at both ends of the 
interval (so called absorbing boundaries) on a lengh of about 450 particles. All these settings are 
summarized in figure [l]. The link with the continuous system is acheived with step h = 1/30, 
for which we shall check that the generated nonlinear structures are continuous objects. 

The expression (0) is inserted in system (D) which is then solved for f n (t) as a system of 
2N coupled first order differential equations in the variable t (N = 15000) with the subroutine 
DOPRI 5. This method allows to launch a plane wave in the system and its efficiency is first 
checked by equating both media to the linear one and we obtain a nice plane wave propagating 
without distortion for very long times (longer than the actual nonlinear experiments) as seen on 
figure §(a) where the parameters of the (homogeneous, linear) medium are c\ = 8, uj\ = 1 and 
the input wave has frequency Q = 1.29 and amplitude A = 1.42. Note in particular that the 
created wave does vanish for large x = nh thus testing the efficiency of the absorbing boundaries 
(the left-hand side absorbing boundary from x = —300 to x = —315 is not represented but 
works as well). 

Second we have verified that, if the right hand side medium is linear, a wave with frequency 
in the phonon gap is indeed totally reflected as seen on figure of |2|(b), in which the amplitude of 

^ee the comment at the end of the paper 
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Figure 1: Grid data for numerical integration 
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Figure 2: Graphs of u(x,t) at t = 250. a) evidence of plane wave generation in a linear 
homogenous medium, b) evidence of total reflection on the stop gap of a linear right- hand- side 
medium. 




the incident wave largely exceeds the threshold amplitude required for nonlinear transmission 
described later. Here the parameters are c\ = 8, u)\ = 1 for the left (linear) medium and 
c\ — 1, u>2 = 3 for the right (still linear) medium. The frequency of the input plane wave is 
Q = 1.29 < 0J2 and its amplitude A = 1.5. 

Results. Thanks to the scale invariance of the system (|1|)(0), we can always fix the parameters 
of one of the two media and vary the others. Herafter we shall then chose the parameters of 
the nonlinear medium to be 

4 = 1, u 2 2 = 3, (8) 

and vary the linear medium (ci and ui) and the incident plane wave Asm(Kx-Qt). We obtain 
the following results. 

1 - An incoming wave with frequency inside the gap of the nonlinear medium (Q < 002) is 
transmitted if its average energy E exceeds a threshold energy E Q . For the choice (|8p we found 
E Q = 1.68. The value of E we refer to is obtained by averaging over one period the hamiltonian 
of the linear chain i.e. 

E = ^A 2 n 2 . (9) 

This energy transmission is accomplished by the generation of kink-antikink pairs and large 
amplitude breathers. The former occur alone when the incident energy equals the threshold 
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Eq. Both breathers and kink- ant ikinks then propagate in the medium. To illustrate this, a 
typical experiment at energy E > Eq is shown on figure ||] where the function u(x,t) is drawn 
at time t = 120 in both media. By looking the picture at different times, we can state that 
the transmitted wave is made of two breathers followed by two kink-antikink pairs, followed 
themselves by a new breather. The various behaviors are then depicted on fig.|] where the wave 

10 | , , , , , , , , 1 
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Figure 3: Solution u(x, t) on the full line att = 120 for a incident wave atVL = 1.29 and energy 
E = 1.8. The parameters of the linear chain are c\ = 8, uif — 1. 

in the nonlinear medium is drawn for different times either for an incident energy E < Eq (total 
reflection), E = Eq (breather generation) and E > Eq (kink-antikink and breather creation). 
Note that the average full width at half maximum of the breather, e.g. the one appearing in 
fig.|] at E — E , contains about 73 particles, hence being a continuous object. 

2 - The found threshold energy E for nonlinear transmission has the following properties: 
it does not depend on the incident frequency Q and it does not depend either on the nature 
of the linear left-hand-side medium (electromagnetic vacuum for io\ = or dielectric medium 
for uj\ 7^ 0). This is illustrated on fig.[| where the dots denote values above which nonlinear 
transmission becomes effective as functions of the incident frequency Vt. 

3 - This transmission process does not occur if the nonlinear system does not possess soli- 
ton solutions. For instance, replacing the nonlinearity sinu by any of its truncated Taylor 
expansions, we do not obtain energy transmission. We illustrate this on fig.|6] where sin u has 
been replaced by its Taylor expansion up to order u 7 . This indicates that this process does not 
result from a shift of the stop gap due to nonlinearity. If such would be the case, one would also 
observe transmission with the Taylor expansion of sin u (and moreover one would also observe 
transmission of radiation as waves, not as nonlinear objects). 

Conclusion By modelizing the scattering of a plane wave incoming onto a nonlinear medium, 
we have discovered the phenomenon that we call nonlinear transmission, and which relies on 
the property of a soliton-bearing system (like sine-Gordon but unlike nonlinear Klein-Gordon) 
of being able to superimpose nonlinearly exponential waves. With respect to gap solitons 
where the incident radiation flows through the medium thanks to the presence of a soliton who 
shifts the gap (or "opens the door"), here the radiation is totally converted into localized large 
amplitude nonlinear objects. These nonlinear objects cannot be described by perturbation 
approaches but require fully nonlinear treatment (e.g. via inverse spectral transform). 

This nonlinear transmission has immediate important implications such as nonlinear ab- 
sorption of radiation (which could play a role in the interpretation of anomalous absorption 
bands), or else nonlinear tunelling, both in continuous and discrete systems. All this issues will 
be displayed in forthcoming works but we want to stress here that the theoretical interpretation 
of the threshold Eq of incident average energy is still to be understood. 
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Figure 4: Representation of nonlinear oscillators u(x,t) for a incident wave at Q — 1.29 < uj 2 
as a function of x for different times. The energy E is varied with the amplitude A of incident 
wave. The parameters of the linear chain are c\ = 8, uf — 1. 



Note that we have checked that this threshold is not due to discreteness effects (Peierls 
barrier to ovecome) by duplicating the numerical simulations with different steps. Varying 
the step size from h = 1/10 to h = 1/100 we have obtained the same threshold of nonlinear 
transmission. 



Comment The nonlinear superposition of exponentials, say exp[£i(a;, t)] and exp[l; 2 (x, t)], is 
obtained by just writing down the two-soliton solution of sine-Gordon under the form 



u(x, t) = 4 arctan 



ai + a 2 e 
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% - a 2 1 + e& + ^ 
where a\ and a 2 are two arbitrary parameters, and where 

2 £j = ( a j + l/ajV 2 /c 2 (£ - Xj) + (a, - l/aj)uo 2 {t - tj) , 



(10) 



'11^ 



with arbitrary constants Xj and tj (note that those values allows us to fix arbitrarily amplitudes 
and phases of both exponentials). It is remarkable that the functions exp[^(x,t)] are precisely 
solutions of the linearized sine-Gordon equation 



c\v xx + u? 2 v 



v a - <*~r rr a;;/' = 
and the formula ( ]TI1| ) indeed furnishes a nonlinear superposition recipe. 



(12) 
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Figure 5: Energy threshold for nonlinear absorption. The parameters of the linear chain are: 
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Figure 6: This figure shows the solution u(x, t) at t = 250 in the case E 
in fig.\Q when sinu is replaced with its Taylor expansion at order u 7 . 



1.89 > E = 1.68 as 



The kink-antikink results from parameters aj real and of same sign (opposite signs produce 
the 2- kink or the 2-antikink solutions). The breather is obtained by assuming complex constants 
a,j with a 2 — a,\. A quite simple exercise allows one to have insight on the creation and 
propagation of such nonlinear structures. In particular the kink and antikink will be seen to 
separate as time runs, as on the numerical experiments in fig.H (use for instance the parameters 



a\ = —0.85, a 2 = —0.9, X\ = 5, x 2 



5, tj 



in the normalized units, i.e. for c 2 = uj 2 — !)■ 



We note finally that the stationary breather furnishes a striking instance of a nonlinear 
superposition of the evanescent waves in (H). Indeed defining 



a\ = a 2 



c 2 n 

-K h I 

0J 2 UJ 2 



(13) 



we obtain for Xj = tj = 



u(x, t) =4 arctan 



n 



l + e 



-2kx 



(14) 



where the relation (f 
stationary breather). 



between k and Q simply means that we have chosen 



1 (for a 



References 



7 



[1] E. Fermi, J.R. Pasta, S.M. Ulam, Studies of nonlinear problems, Los Alamos Sci. Lab. 
Rep. (1955) LA-1940, reproduced in: A.C Newell (Ed.) Nonlinear wave motion, AMS 
Lect Appl Math 15 (1974) 

[2] N.J. Zabusky, M.D. Kruskal, Phys Rev Lett 15 (1965) 240 

[3] C.S. Gardner, J.M. Greene, M.D. Kruskal, R.M. Miura, Phys Rev Lett 19 (1967) 1095. 
V.E. Zakharov, A.B. Shabat, Sov Phys JETP 37 (1972) 62. M.J. Ablowitz, D.J. Kaup, 
A.C. Newell, H. Segur, Stud Appl Math 53 (1974) 249. 

[4] G.L. Lamb, Jr., Elements of Soliton Theory, J. Wiley and Sons (New York 1980). M. 
J. Ablowitz, H. Segur, Solitons and the Inverse Scattering Transform, SIAM (Philadel- 
phia 1981). F. Calogero, A. Degasperis Spectral Transform and Solitons, North Holland 
(Amsterdam 1982). R.K. Dodd, J.C. Eilbeck, J.D. Gibbon, H.C. Morris, Solitons and 
Nonlinear Waves, Academic Press (London 1982). M. J. Ablowitz, P. Clarkson, Solitons, 
Nonlinear Evolutions and Inverse Scattering, Cambridge Univ. Press (Cambridge 1992). 

[5] M. Remoissenet, Waves Called Solitons, Springer (Berlin 1999). A.C. Scott, Nonlinear 
science, Oxford Univerity Press (New York 1999). 

[6] S.A. Gredeskul, Yu S. Kivshar, A. Sanchez, L. Vazquez, Phys Rev A 45 (1992) 8867. R. 
Scharf, Yu S. Kivshar, A. Sanchez, A.R. Bishop, Phys Rev A 45 (1992) 5369. G. Reinisch, 
Int J Mod Phys B 6 (1992) 3395. Yu S. Kivshar, L. Vazquez, Phys Rev Lett 77 (1996) 
582. C.R. Willis, Phys Rev E 55 (1997) 6097. A. Sanchez, A.R. Bishop, SIAM Rev 40 
(1998) 579. L. Freidland, Phys Rev E 57 (1998) 3494. G. Cohen, Soliton interaction with 
an external travelling wave, preprint flhttp: / / arXiv.org/ ps / patt-sol/ 9903006|) . 



[7] M. Peyrard, M.D. Kruskal, Physica D 14 (1984) 88. R.S. MacKay, S. Aubry, Nonlinearity 
7 (1994) 1623. D. Cai, A.R. Bishop, N. Gronbech- Jensen, Phys Rev Lett 72 (1994) 591. 
S. Aubry, Physica D 103 (1997) 201. S. Flach, C.R. Willis, Phys. Rep. 295 (1998) 182. 
D. Hennig, G.P. Tsironis, Phys Rep 307 (1999) 333. T. Dauxois and M. Peyrard, Phys 
Rev Lett 70 (1993) 3935. M. Peyrard, Physica D 119 (1998) 184. 

[8] S.L. McCaU, E.L. Hahn, Phys Rev 183 (1969) 457 . J.C. Eilbeck, J.D. Gibbon, P.J. 
Caudrey, R.K. Bullough, J Phys A 6 (1973) 1337. G.L. Lamb Jr, Phys Rev A 8 (1974) 
422. M.J. Ablowitz, D.J. Kaup, A.C. Newell, J. Math. Phys. 15 (1974) 1852. D.J. Kaup, 
Phys Rev A 16 (1977) 704. 

[9] F.Y.F. Chu, A.C. Scott, Phys. Rev. A 12 (1975) 2060. K. Druhl, R.G. Wenzel, J.L. 
Carlsten, Phys Rev Lett, 51 (1983) 1171. D.C. MacPherson, R.C. Swanson, J.L. Carlsten, 
Phys Rev A 40 (1989) 6745. J.C. Englund, CM. Bowden, Phys Rev A 42 (1990) 2870, 
and Phys Rev A 46 (1992) 578. C. Claude, J. Leon, Phys Rev Lett 74 (1995) 3479. C. 
Claude, F. Ginovart, J. Leon, Phys Rev A 52 (1995) 767. 

[10] J. Leon, A.V. Mikhailov, Phys Lett A 253 (1999) 33. M. Boiti, J-G. Caputo, J. Leon, F. 
Pempinelli, Inverse Problems, 16 (2000) 303 

[11] H. Steudel, Physica 6D (1983) 155. D.J. Kaup, Physica 6D (1983) 143, and Physica 19D 
(1986) 125. 

[12] G.P. Agrawal, Nonlinear fiber optics, Academic Press (London 1989) 
[13] D.L. Mills, Nonlinear Optics, Springer Verlag (Berlin 1991) 



R 



[14] I.R. Gabitov, V.E. Zakharov, A.V. Mikhailov, Theor Math Phys 63 (1985) 328, and Sov. 
Phys. JETP, 59 (1984) 703. See also J. Leon, Phys Lett A 170 (1992) 283; Phys Rev 
A47 (1993) 3264; J Math Phys 35 (1994) 3054 

[15] A. Hasegawa, F. Tappert, Appl Phys Lett 23 (1973) 142. L.F. Mollenauer, R.H. Stolen, 
J.P Gordon, Phys Rev Lett 45 (1980) 1095 

[16] W. Chen, D.L. Mills, Phys Rev Lett 58 (1987) 160. D.L. Mills, S.E. Trulliger, Phys Rev 
B 36 (1987) 947. W. Chen, D.L. Mills, Phys Rev B 36 (1987) 6269. 

[17] E. Lidorikis, Q. Li, CM. Soukoulis, Phys Rev B 54 (1996) 10249. 

[18] C. Martjin de Sterke, J.E. Sipe, Phys Rev A 38 (1988) 5149. 

[19] Photonic band gaps and localization, CM. Soukoulis (Editor) NATO ARW, Plenum (Ney 
York 1993). Photonic band gap materials, CM. Soukoulis (Editor) NATO ASI, Klumer 
(Dordrecht 1996). Photonic crystals: molding the flow of light, J.D. Joannopoulos, R.D. 
Meade, J.N. Winn, Princeton Univ. Press (Princeton 1995). 

[20] F-X. Hugot, J. Leon, Inverse Problems 15 (1999) 701 

[21] J-G. Caputo, A.C Newell, M. Shelley, in Integrable Systems and Apllications, Eds. M. 
Balabane, P. Lochak, D.W. McLaughlin and C. Sulem, Springer (Berlin 1990). 



Q 



